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Abstract. The planar elliptic extension of the Laplacian growth is, after a 
proper parametrization, given in a form of a solution to the equation for area- 
preserving difFeomorphisms. The infinite set of conservation laws associated 
with such elliptic growth is interpreted in terms of potential theory, and the 
relations between two major forms of the elliptic growth are analyzed. The 
constants of integration for closed form solutions are identified as the singu- 
larities of the Schwarz function, which are located both inside and outside the 
moving contour. Well-posedness of the recovery of the elliptic operator gov- 
erning the process from the continuum of interfaces parametrized by time is 
addressed and two examples of exact solutions of elliptic growth are presented. 



1. Introduction 

Several moving boundary processes, such as solidification [1], electrodeposition 
[2], viscous fingering p], and bacterial growth to name a few, can be reduced, 
after some idealizations, to the Laplacian growth, which can be described as follows: 

(1.1) V{0 = dnGD(t){^,a). 

Here V is the normal component of the velocity of the boundary dD{t) of the 
moving domain D{t) C M'*, ^ G dD{t), t is time, 9„ is the normal component of the 
gradient, and GD(t){£,, o) is the Green function of the domain D{t) for the Laplace 
operator with a unit source located at the point a € D{t). 

In two dimensions this equation can be rewritten as the area-preserving diffeo- 
morphism identity 

(1.2) Q{ztz^) = l, 

where z{t, </>) := dD{t) is the moving boundary parameterized by ^ G [0, 2it\ and 
conformal when analytically extended in the region SJc/) < (Sj [6] . The equation 
()1.2p possesses many remarkable properties, among which, the most noticeable ones 
are the existence of an infinite set of conservation laws: 

(1.3) C„ = / z'^dxdy, 

JD{t) 
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where n runs over all non- negative j7j (non-positive [8j) integers in the case of a 
finite (infinite) domain D(t), and an impressive list of exact time-dependent closed 
form solutions [2l\. For a beautiful interpretation of conserved quantities C„ as 
coefficients of the multi-pole expansion of the fictitious Newtonian potential created 
by matter uniformly occupying the domain D{t) see, e.g., |21) . 

It was established in [22. that the interface dynamics described by (|1.2p is equiv- 
alent to the dispersionless integrable 2D Toda hierarchy [25 j constrained by the 
string equation. Remarkably, this hierarchy, being one of the richest existing inte- 
grable structures, describes an existing theory of 2D quantum gravity (see the com- 
prehensive review [23] and references therein). The work [22] generated a splash 
of activity in apparently different mathematical and physical directions revealing 
profound connections between Laplacian growth and random matrices j24l ■ the 
Whitham theory [25j . and quadrature domains [26j . 

In this paper we present a natural extension of the Laplacian growth, where the 
Green function of D{t) for the Laplace operator in the RHS of (1.1) is replaced 
by the Green function of a linear elliptic operator, 

(1.4) L = V • (A(x)V) -u(x), A(x)>0, x e M"*. 

Such a process, which is natural to be named an elliptic growth, is clearly much 
more common in physics than the Laplacian growth. 

Consider, for instance, viscous fingering between viscous and inviscid fluids in 
the porous media governed by Darcy's law 

(1.5) v = -AVp, 

where A is the filtration coefficient of the media and p is the pressure (equal to the 
Green function, G]j(t), defined in (1.1) in most of the cases of interest for us). One 
can easily imagine a non-homogeneous media where the filtration coefficient A is 
space-dependent. Such examples of elliptic growth, where the elliptic operator L 
has the form of the Laplace-Beltrami operator, L = V • AV, and A is a prescribed 
function of x, will be called an elliptic growth of the Beltrami type. It is clear that all 
moving boundary problems other than viscous fingering with a non-homogeneous 
kinetic coefficient A fall into this category. 

From a mathematical point of view this process is the Laplacian growth occurring 
on curved surfaces instead of the Euclidean plane. In this case the Laplace equation 
is naturally replaced by the Laplace-Beltrami equation, and A (that can be a matrix 
instead of a scalar as it is in our case) is related to the metric tensor. There are 
several works addressing the Hele-Shaw problem on curved surfaces and we will 
mention below those few related to the integrable mathematical structure of elliptic 
growth. 

Another major source of examples of elliptic growth is related to screening ef- 
fects, when u ^ 0, while A is constant in (|1.4p . The simplest example of this kind 
is an electrodeposition, where the field p is the electrostatic potential of the elec- 
trolyte. It is known that in reality electrolytes ions are always locally surrounded by 
a cloud of oppositely charged ions. This screening modifies the Laplace equation for 
the electrostatic potential by adding to the Laplace operator the negative screening 
term, —u{x), which stands for the inverse square of the radius of the Debye-Hukkel 
screening in the classical plasma [27]. For the homogeneous screening m is a (pos- 
itive) constant, so the operator L becomes the Helmholtz operator, while for the 
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non-homogeneous case, when u is not a constant, L is a standard Schrodingcr op- 
erator. Motivated by this example, we wiU call the moving boundary problem for 
L = — u an elliptic growth of Schrodinger type. 

We show that these rather general types of elliptic growth still retain remarkable 
mathematical properties, similar to those possessed by the Laplacian growth. A 
mixed case with a non-constant A and non-zero u also shares similar properties but 
is less representative in physics and can always be reduced to one of the two former 
types of elliptic growth by a simple transformation described later on in the article. 
For completeness we shall indicate another class of elliptic growth when the fluid 
density p changes in space while it is constant in time. This happens for instance 
when porosity (fraction of porous media accessible for fluid) is space-dependent. In 
this case the continuity equation for incompressible fluid in porous media has the 
form 

(1.6) V(p(x)A(x)Vp) = 0, 

while p.5p still holds. This case presents an additional extension of the elliptic 
growth related to potential theory with a non-uniform density, as will be shown 
below. 

In prior works on elliptic growth an infinite number of conservation laws, re- 
garded as extensions of (|1.3p . were identified in [211 [28]. Also an integrable exam- 
ple in 2D, which corresponds to a very special choice of the conductivity function, 
A(x), was explicitly constructed in [29]-[3T]. The elliptic growth in these works was 
reduced to the well-known Calogero-Moser integrable system. 

The present article contains several new results on elliptic growth and reviews 
the known conservation laws from a slightly novel perspective. It is organized as 
follows: 

- In Section 2 we interpret (|1.2p as the equation of the area-preserved diffeomor- 
phism in 2D and analyze its connections with the Laplacian growth. 

- Section 3 contains the definition of elliptic growth, the conservation laws for 
this process and recasts the latter in terms of the inverse non-Newtonian potential 
theory. 

- In Section 4 we obtain the equation (|1.2p for elliptic growth of the Beltrami 
type by introducing the function q, conjugate with respect to p defined in (|1.4p . 
which plays the role of a stream function for the incompressible fluid; furthermore, 
we obtain the Beltrami equation for the function p -f iq. 

- Section 5 is devoted to analyzing connections between the elliptic growth of 
the Beltrami and Schrodinger types, elucidating the difficulties of parametrization 
of the interface for the Schrodinger type. 

- In Section 6 we reformulate the elliptic growth in terms of the Schwarz function 
of the moving interface. 

- Section 7 addresses a well-posedness of a recovery problem for the operator L 
from the continuum of moving interfaces parameterized by time. 

- In Section 8 we discuss Herglotz' theorem as the main device to generate 
exact solutions and present two examples of the exact closed form solutions of the 
elliptic growth. We also identify the constants of motion of these solutions as the 
singularities of the Schwarz function of the moving contour. 

- Section 9 contains brief conclusions. 
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Due to the fact that one of us is labelled as a theoretical physicist and following 
the customs of the physics community the references do not appear in alphabetical 
order. 
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1.1. List of notations and conventions. We collect below a few basic definitions 
and notations used throughout the text. 

V2 = A, V/ = grad/; 

V(U) = V • U = divU, where U is a vector field; 

denotes the class of real analytic functions; 
h - ^■ 

n stands for the outer unit normal to the moving boundary T — T{t)\ 
£ denotes the arc length on the boundary F; 
dA=^dArea= ^^^; 

7r 7r ' 

an analytic Jordan curve means a smooth Jordan curve which admits a real 
analytic parametrization. 

2. Area preserving diffeomorphisms 

This section contains some immediate implications of the equation of area pre- 
serving diffeomorphisms related to the parametrization of an analytic Jordan curve. 
Later on we shall see that this, apparently innocent, Jacobian identity plays an im- 
portant role in the study of moving boundaries governed by elliptic growth. 

2.1. Fourier expansion. Consider the equation 
(2.1) ^{zlz,) = 1, 

where t e [0, T] is a non- negative variable (usually identified with time), while 
q € [0, 27r] is the parameter along the Jordan analytic curve Ct- Equation (|2.ip 
can be interpreted as an area preserving property: that is the Jacobian of the 
transformation 

{t,q) 1-^ {x,y), where z{t,q) = x + iy, 

is equal to 1. 

To be more precise, for a fixed t, we assume that the 27r-periodic real analytic 
map 

2(t, •) : [0,27r] — > C 

is an embedding, and its range is denoted by Ct- We denote by D{t) the interior 
of the Jordan curve Ct- 

In view of the smoothness hypothesis imposed on z{t,q) we can expand the 
function z{t, q) in a Fourier series 



(2.2) z{t,q)^Y,''k{t)e' 



kq 
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We will assume that the dependence t i— > z{t,q) is C^. Also, we put w = e*'', so 
that we can rewrite (2.2) as 

oo 

z{t,w) = y^flfciQw''. 

— oo 

By the analyticity assumption, there exists e,0 < e < 1, so that the above 
Laurent series is convergent in the annulus 1 — e< < 1 + e. 

Due to the real analyticity of the map z, the Fourier series for z and its derivatives 
are absolutely and uniformly convergent, whence 

oo 

— OO 

and 

oo 

Thus, equation (|2.ip becomes 



-oo 



1 = ^{ztZq) = i '^{a.nkak + dknan)e 

k,n 



2 " 



OO ^ OO 

m — — 00 n— — 00 

By equating the coefficients we find 

^ 00 

n— — 00 

and 

00 

0=- ^ {{n + m)anan+m + nan+mO^) J 



2 

n— — 00 



whenever m ^ 0. 
Since 

Area(i.(t)) = / ^ = / W:^)^^'^ = /"iM)'^^^'^^ 



iCt 2« J\w\=i ' dw 2i Jo ' dq 2i' 
we derive the following remarkable identity. 



Proposition 2.1. Under the assumption i2.1\) . the family of domains bounded by 
the curves z(t.,q), < q < 2it, satisfy 

dATe&{D{t)) 



dt 



= 2tt. 



By regarding z now as a function of t and the complex variable w = e''^, we 
obtain, along the curve Ct, the following: 

dz dz dw 
oq aw oq 

Denote, for the sake of simplicity, z' = Zw Then, the master equation (|2.ip 
becomes 

3?(ii;z'i) = 1. 
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Since z{t, •) is analytic in the annulus < 1 — e < 1^1 < 1 + e, the foUowing result 
follows. 

Proposition 2.2. Under the above assumptions 

wz'{t,w)i\t,l/w) + -z'\t,-)z{t,w) =2, 1 - e < < 1 + e. 
w w 

(As usual, for a complex analytic function h{w), we denote by 

h^{w) = h{W), 
obtained from h by conjugating its Taylor coefficients.) 

2.2. Analytic parametrization. The most studied case of the Laplacian growth 
process requires an additional analyticity assumption. We devote the present sub- 
section to this scenario. Assume that for all t E [0, T] the negative Fourier coeffi- 
cients vanish, i.e., 

(2.3) ak{t) = 0, fc < 0. 

It is not difficult to see that this will hold the whole evolution, t € [0,T]. The 
equation (|2.3p simply means that for a fixed t the function z{t,w) extends analyti- 
cally to the unit disk w g D. Since z(t, •) is a homeomorphism from the boundary 
T = to the curve dD{t), the argument principle implies that 

z{t,-) : D — > D{t) 

is a conformal mapping. Let w ~ 'i>{t,z) denote the inverse conformal mapping. 
Since a linear transformation z az+(3, |a| = 1, leaves the equation (|2.ip invariant, 
we can assume without loss of generality that z{t, 0) = and that p{t) = z'{t, 0) > 0. 
To distinguish this case from the general case considered in the previous subsection 
we shall denote dD{t) by V{t). The function p{t,z) = log |vl'(t, z)| is, up to a 
constant factor, the Green function of the domain D{t), with the source at z = 0. 
This means that p{t, .) is the unique harmonic function in the punctured domain 
D{t) \ {0} having zero boundary values on r{t) and such that p{t,z) — log|z| is 
harmonic at z = 0. 

Moreover, the harmonic conjugate function arg^'(t, z) is, up to an additive con- 
stant, equal to q{z), z e T{t). 

In other words, for a fixed value of the parameter t, we have: 

V^p{t,-)^27rS{-), in D{t), 
Pit,-)\m=0: 

and 

{Qvit, ■),-qxit, •)) = {Px{t, ■),Pyit, •))■ 

Consequently, the normal velocity of the boundary equals 

y^ dq{t,-) ^ dp{t,-) 
d£ dn 
So, by the area conservation property, we have 

Vdz Ad£ — dz A dq. 

These equations define a specific dynamics of planar boundaries known as Lapla- 
cian growth. For recent guides to the mathematics and physics behind Laplacian 
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growth we refer to the volume [32] and the survey l26j . We will return to this case 
after discussing the geometry of the moving boundaries. 

2.3. The Schwarz function. An important tool for studying the changing geom- 
etry of the moving boundaries is the Schwarz function [ 33] , [45] . Up to the complex 
conjugation it is simply the (local) Schwarz reflection with respect to an analytic 
curve. 

On the real analytic smooth boundary T{t) of D{t) we introduce the Schwarz 
function 

z = S{t,z), 

where S is analytic in the variable z. The domain of definition for S(t, .) is at 
least a tubular neighborhood of T{t), although the function may possess analytic 
extensions to much larger sets. For instance, the Schwarz function of a disk centered 
at z = a and of radius r is the rational function 

r2 

S(z) = a+ . 

z — a 

If a polynomial P{z,'z) vanishes on T{t), then, necessarily, the associated Schwarz 
function satisfies the algebraic equation 

P{z,S{t,z)) = 0, zeT{t). 

Proposition 2.3. 34J The normal velocity of the boundary satisfies 

V 



2iVS; 

with the proper choice of the branch of the square root, so that \/\fSl = dz/dl 
along Tit). 

Proof. By taking derivatives with respect to t we have 

zT = 5t + SzZt. 
When restricted to the boundary curve, 

S ^ — 

^ dz 

is a complex number of modulus one. 

Fix a single- valued branch of \fSl along dD{t). This is always possible since 
^1 \fSl equals to the unit tangent vector to dD{t) and, hence, is single valued near 
dD{t). Then the above equation becomes 

— St /TT Zt 



or, equivalently. 



St „ Zt 



/s: 



Since 



is the unit tangent vector, then 



1 dz 



St 



= 2iV 

□ 
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Theorem 2.4. There exists a multivalued analytic function W{t,z), defined in a 
neighborhood of dD(t), with the property 

St = dM, 

and such that is constant along dD{t). 
Proof. From the above computations we find 

-2i 



St = -2i^S,V = -j=V, 

\ Sz 

whence the vector St is collinear with the complex conjugate of the normal to T{t). 
Moreover, rewriting the last equation in the form 

^t^dW 

we infer 

\St\n = 'St:^dW = ^{m^), 
where n is the normal to dD{t). Hence, it follows that the boundary, dD{t) is a 
level set of IkW. □ 

2.4. Laplacian growth. In this subsection we merely illustrate few classical ob- 
servations related to the consequences of the dynamics (|2.ip under the analyticity 
assumption ak{t) = 0, A: < 0, i € [0,T]. That is, we assume again that the 
parametrization •) of the curve V{t) analytically extends to the interior of the 
unit disk B and will use intensively the Schwarz function techniques. 

By returning to the notations introduced in Section 12.21 we can identify the 
complex potential W with the multivalued function C,{t^ z) = p{t, z) + iq{t^ z) — 
log '^{t, z) and then study the analytic extension of the Schwarz hmction S{t, z). 

Theorem 2.5. [34 For every t, there exists a tubular neighborhood U ofT{t), such 
that 

(2.4) St{t,z) = 2Ut,z), z&U. 

Note that the function C, is multivalued and analytic in the punctured domain 
D{t) \ {0}. Its derivative is therefore meromorphic there with a simple pole at 
z — Q and the residue equal to 1. 

Proof. Indeed, according to Proposition 12. 31 we have along T{t): 

St ^ dp 
2iVK On 

so 

^ _ 2i dp _ 2i{dp/dn)d£ _ 
* dz/ di dn dz 



^ {idq/d£)d£ _ ^ d{p + ig) / dl)d£ _ ^dC 
dz dz dz 



□ 



As simple as it looks, equation (|2.4p has surprising consequences. In order to 
unveil them, we start with the known Plemelj-Privalov-Sokhotsky formula applied 
for a fixed t to the function S{t, •). Let 

S±(t,z)^— / zeD(t), respectively z £ C\D(t). 

2711 Jp(t) a~ z 
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Then 

(2.5) S{t,z) ^ S+{t,z) - S-{t,z), zeU, 

where U denotes, as before, a neighborhood of T(t). Similarly we decompose the 
fmiction (z{t,z) and find 

{C)_it,z) = —. 

z 

From the uniqueness of the above decompositions and (|2.4p . we infer 
or, for z ^ D{t), 



S-{t,z)t^^, 



d 1 r ada _-2 
dt 2m Jr{t) f ~ ^ ^ ' 

that is 

d r dA{a) _ -2 
dt Joit) <^ - z z ' 
By integrating against a polynomial / along the circle \z\ = i?, with R sufficiently 
large, we find the following general identity. 

Proposition 2.6. // the parametrization z{t,w) of the boundary of the domain 
D(t) extends analytically to the interior of the unit disk D := {w : \w\ < 1}, then 
for every polynomial f{z) the following identity holds: 

^ I f{z)dA{z) = 2/(0). 

«t JD{t) 

Equivalently, S-{t,z) — + h{z), where h{z) is an analytic function in the 
neighborhood of C \ D{t), vanishing at infinity, and independent of t. A simple 
application of Cauchy's formula now yields that there exists a complex valued 
measure /i supported on a compact set K C D{t), independent of t (as proved 
above) and such that 

'2t 1 f d^{a) 



S-it,w)^ / zeC\D{t). 

z n a - z 

By repeating the above calculations wc find 

dA{a) _-2t If d^{a) 



ID(t) - ^ ^ 1^ JK <J - ^ 

and, consequently, the following quadrature identity follows. 

Corollary 2.7. Under the same hypotheses as in Proposition 2.6, for every poly- 
nomial f we have 

(2.6) / f{z)dA{z) = 2tf{0) + - [ f{z)dii{z). 

JD(t) Jk 

Simple examples show that the measure /x is not unique. If one insists that the 
supporting set K is "minimal", and the representing measure jjL is positive, then 
one can prove in most interesting cases the uniqueness of fx. The case of quadrature 
domains -D(t), corresponding by definition to a positive finite atomic measure /z, 
is by far the best understood from the constructive point of view. In this case 
the conformal mappings z(t, w) are rational. Examples, a discussion of the alluded 
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uniqueness and further details and references can be found in the coUection of 
articles [35], cf. also [42], citeShapiro. 

3. Elliptic growth 

Guided by Laplacian growth as a prototype, we introduce in this section the 
elliptic growth phenomenon mentioned in the Introduction. It is surprising to see 
that many features of Laplacian growth persist and yet sharp differences occur. Let 
us start the formulation in arbitrary dimension d for a possibly multiply connected 
domain D{t) in R'' with many sources, but later on we will focus on a homotopically 
trivial 2D case with a single source at the origin in more detail. 

Consider a family D{t) of bounded domains in with smooth analytic bound- 
aries. Moreover, dependence of D{t) on t is assumed to be real analytic (in the 
sense of a chosen parametrization) as well. 

Let G be an open set containing as relative compact subsets all D{t), —1 < t < 1, 
and let A : G — > (0, oo) be a real analytic function. We consider the elliptic (non- 
positive) differential operator 

L = V • AV = div(Agrad). 

As we noted in Section 1.1, when there is no danger for confusion we shall omit the 
dot in the notation, and write, for example, A = V^. 

The moving boundary problem with N sources Sk at e D{t) is the following: 

Given -D(O), find domains D{t) satisfying the system of equations: 

N 

Lp = ^Sfe(5(x-Xfe) in D{t), 
fc=i 

P\dDit) = 0, 

V = XdnP on dD{t). 

Note that the first two conditions simply assert that p is the linear combination 
of the Green functions for the operator L of the domain D{t)^ with singularities at 
Xk, while the third condition determines the dynamics of the moving boundary. 

Theorem 3.1. For every function £ C'^(G) satisfying Lxp ~ 0, we have 

d r ^ 

Here, dVol stands for Lebesgue measure on D{t). 

Since the constant function -0 = 1 is annihilated by the operator L, the above 
formula implies 

d ^ 
-Voli?(i)=^s,, \t\<l. 

fe=i 

Thus, in this moving boundary process, the volume is still proportional to time. 
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Proof. Let dT denote the surface element on each connected component of the 
boundary of D{t). Then, wc have: 

d_ 
di 



D{t) JaD(t) JdD(t) 



I [ipWIp - pXS/ip) -ndV^ I V • (V'AVp - pAVV')dVol 

JdD(t) JdH) 

f ^ 

/ [i^Lp - pL4']dVo\ = Sfc'(/'(xfc). 
JdU) 



D(t) k=l 

□ 

Corollary 3.2. In the case when the domains D{t) are all homeomorphic to a ball 
and contain a single source Si > 0, the moments 



C{^) = J VdVol, = 0, 

determine the domains D(t) (locally int). 

Proof. Indeed, it is sufficient to consider a single moment C(l). As remarked 
earlier, 

^^^=/ VdT^s,>0 

dt JdD(t) 

and the corollary follows, after observing that the family D{t) is increasing with 
respect to the ordering by inclusion. □ 

Remark. Theorem 13 . II and the Corollary extend word for word to more general 
elliptic operators 

L = div (A grad) — w(x), 

where the matrix A = (A^j (x)i.j=i,2) is uniformly elliptic on the domain G and all 
the coefficients A^.j, m(x) are assumed to be real analytic in G and u > 0. The only 
modification needed is that in the last boundary condition in (1.1) where one ought 
to require V — A( Vp) -n. The existence of the Green function p for such operators is 
well known [36l[37]. The fact that > on dD{t) then follows from the maximum 
principle and Hopf's lemma which hold for such operators cf. [551 [571 

Assuming that the sources strengths Skit) depend on time we then obtain another 
notable corollary of Theorem l3.1l The functionals i^d^ do not depend on Sk{t), 

but only on the value of the integral s(t) dt [Ml HD US- 

The functionals /o^j-) ^dYo\ have a remarkable potential theoretic interpretation 
[211 |28j . Indeed, imagine that a domain D{t) is occupied by matter with a unit 
density, which creates the potential $, governed by the Poisson's equation 

= XD(t), 

where xd is the characteristic function of the domain D. A solution of the last 
equation is 



<I>(x)= / Go(x,y)dVol(y), 
Jdh) 



12 



KHAVINSON, MINEEV, AND PUTINAR 



where Go(x, y) is the fundamental solution for the operator L. In important par- 
ticular cases which are relevant for physical applications (for instance, for the 
Helmholtz operator A — 1) it is possible to expand G'o(x, y) into the series 

Go(x,y) = ^V^„(x)^„(y), x i D{t),y e D{t), 

n 

where {"iAn} and {"iAn} are bases of the null space of L in D{t) and its complement 
in R'^ respectively. Then, assuming commutativity of summation and integration, 
we obtain 

<i>(x) = V^„(x) / ^„(y)c?Vol(y). 
JD{t) 

Therefore, we have obtained the functionals introduced in the Theorem 13.11 as the 
coefficients of the multi-pole expansion of the non-Newtonian potential given in 
a far field. We would like to add that the gradient of $ is a generalization of 
the Cauchy transform for the domain D{t) - the notion that was so useful in the 
Laplacian growth and the related field of quadrature domains. 

In the case of elliptic growth with nonhomogeneous density mentioned in the 
Introduction, one should modify the formulation given in the beginning of this 
section by adding a positive space-dependent factor, p, namely, 

A -> pA; 

under these assumptions we still encounter an infinite set of conservation laws 
similar to the previous case, when p = 1, namely: 

— I -dVol = V9fcV'(xfc). 
diJoit) P ^ 

The elementary proof is not included here. From the point of view of potential 
theory this corresponds to the case of occupation of the domain D{t) by non- uniform 
matter with density l/p(x). Some aspects of this case in 2D were discussed in [39], 
also cf. [13]. It is clear that the inverse potential problem of recovery of D{t) is 
considerably more difficult in this situation. 

4. The conjugate function 

From now on we return to a planar case with a single source of strength si = 2tt 
and assume that all domains D(t), < t < T, are simply connected. Define, 
using the notations from the previous section, a (multivalued) conjugate function 
q G C^{D{t)) by 

(4.1) qy = Xpx 

(4.2) = -Xpy. 

Accordingly, 

V • ivq = 0, in D{t). 
A 

In most computations below t is fixed. However, we stress that by its very definition, 
the function q depends on t also: q{z,z) = q{t; z,z). We hope that omitting t in the 
notations of q will not confuse the reader. 
We have 

V = XdnP = q. 
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Let z{t,£) be the parametrization of the contour T{t) ~ dD{t) by an arc-length 
Then a right angle rotation of the unit tangent vector gives 

n — —izg, 

hence the normal component of the boundary velocity is 



V ^ Zf {~ize) = ^{zt{-izi)) = "^{ztze) = deq, 

where, as before, subscripts stand for partial derivatives. 

Since A is positive, then deq > along r(i), therefore q{t, .) can equally well 
parameterize the boundary dD{t) = V{t), in which case we can rewrite the above 
equation as 

3(ztZq) = 1, on dD(t). 



Lemma 4.1. The variation of q along the curve T{t) is equal to 27r. 
Proof. Indeed 

var q\r{t) 



iT{t) de Jr{t) dn 

/ AVp -nde^ / W{X'Vp)dArea = 2tt. 
Jr(t) JDU) 



-(f) JD{t) 

Introducing the multivalued function 

C,{z,z) ^p + iq, 

we have 

iXdp — dq, 

or, still, 

(4.3) Be - ^9C, 

which is a form of the Beltrami equation 44J. In terms of a new variable 

uj = \/Xp + i-^, 

Va 

the equation (|4.3p takes the canonical Carleman form 

(4.4) Boj = {d\ogVX)LJ. 
An implicit solution of this equation is found from 



□ 



( , f aiog(yA(c^ ^(c^ 

uj{z,z) = F{z)exp — ^dA(C), 

where F{z) is analytic in D{t) (cf. |44j for more details). 

Thus the moving boundary problem of finding D(t) can be reformulated as a 
Dirichlet boundary value problem: 

Given the weight X, find a function <^(z,z) (or, uj{z,z)) satisfying the Beltrami 
(or, Carleman) equation above and subject to the boundary condition SR^ = (or, 
^uj = 0). 
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5. Elliptic growth of Schrodinger type 

The previous section was devoted to the elhptic growth of the Behrami type. 
In this section, we will consider the elliptic growth of Schrodinger type, which, as 
already mentioned in the introduction, is related to the theory of the Schrodinger 
operator. To be specific, we consider the problem: 

Given a domain -D(O) find domains D{t) , satisfying the system of equations: 

N 

LP = {V^-u)P = ^qk6{x-Xk) in D{t), 
fc=o 

P = on dD{t), 

V = dnP ondD{t). 

As has already been demonstrated, this problem has an infinite set of conservation 
laws which are time derivatives of integrals of null vectors of the operator L. These 
integrals have the potential-theoretic interpretation discussed in Section 3. Let us 
pose the following question. Does there exist in this case a function Q "conjugate" 
w.r.t. to P in the sense that P and Q are connected via some generalized Cauchy- 
Riemann equations? And, if such Q exists, can it be used as a parametrization of 
the moving contour dD{t) similarly to how it was used in the case of the elliptic 
growth of the Beltrami type in the previous section? 

The answer to the first question is 'yes', and to the second one - 'no'. One can 
see this from the generalized Cauchy-Riemann conditions which connect p and q in 
the Beltrami case. Namely, 

Xd^p = dyQ, 

Xdyp = -d^q. 
These formulae suggest the substitution 

P 

q = VXQ. 

Thus, the new functions P and Q are connected via the system of linear euqa- 
tions: 

d^P - P9x(log VA) = dyQ + dy{\og VA), 
dyP - P9y(log VA) = -d^Q - a,(log VA). 

Differentiating the first equation w.r.t. to x, the second one w.r.t. to y, and then 
adding them, one obtains 

(V=^ - u)P = 0, 
(V^ - v)Q = 0, 

where 

_ V"(Ai/^) 

V2(A-V2) 

i; = 

A-V2 

This simple transformation known as the "removal of the first derivative from linear 
differential equations of the second order" and, also, closely related to supersym- 
metry in physics, joints together the two major types of elliptic growth. 
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Now let US show that the function Q, unlike the function q cannot, in general, 
provide a parametrization of the contour. Indeed, it was shown in the previous 
section that q can serve as a parametrization since it is a monotonically increasing 
function of the arc-length along the contour. Since q it is now clear that 

Q, generally speaking, is not monotone along the interface because of the space 
dependent factor A~^/^. 

As one can easily see, a/A solves the same Schrodinger equation as P, namely 

(V^ - w)\/A = 0, 

while the function l/\/A solves the same Schrodinger equation as Q, namely 

6. An inverse problem 
We address below the following natural question: 

Is it possible to have the same "movie" {t,T{t)), t G [0,T] governed by the ellip- 
tic growth dynamics with different weights A ? 

By studying a particular example we shall demonstrate that, indeed, such 
non-uniqueness may take place. To fix the ideas, assume that the elliptic growth 
dynamics, as specified above has the property that the conformal map z{t, •) onto 
a neighborhood U of T{t): 

z{t, •) : {w; 1 - e < < 1 + e} — >U 

extends analytically to the unit disk \w\ < 1. Then, as we saw earlier, the equation 
(|2.H) once again governs the whole evolution of T{t) but this time the evolution of 
r(t) is the Laplacian growth. In particular, the parameter along each boundary 
satisfies q{t,z) = Qlogw{t, z), that is V^g = at all points except the isolated 
singularity. On the other hand we have started with the assumption VA^^Vg = 0. 
Hence 

A~^V\-A"2(VA)-(Vg) = 0, 

so 

(VA) • (Vq) = 0. 
The function p{t, z) also, by its very definition, satisfies 

(Vp) • {Vq) = 0. 

Since Vp never vanishes, at least in a neighborhood of r(i), a functional dependence 

\{z)^ f{t,p(t,z)) 

must occur. This simple computation gives rise to the following counterexample. 

Let us consider the simplest Laplacian growth process of expanding concentric 
discs 

D{t) = D{Q,V2i), t>0, 

with associated functions 

z) = log— =, g(<,z) = argz. 
V 2t 
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Let r = \z\ and choose any positive, smooth radial function A(r), e.g., A(r) = expr. 
Then, everywhere in C \ {0}, we find 

VA(r)-^Vg = \{r)-^Aq - A'^VA • = 0. 

Now we can solve the problem with a rotationally symmetric function p{t,r): 

V\{r)Vp = 27r5(0), p{t, \/2t) = 0, 

so that p, q are L-conjugate, where L = VAV. In particular, the normal velocity of 
the boundary is given by 

dq dp 

thus it is independent of the choice of A. That is, the same movie can be described 
as the Laplacian growth with A = 1 and as elliptic growth with any other positive 
weight which is rotationally invariant. 

Note that a point on the boundary r{t) = dD{0, is parametrized as 

z{t,q) = V2ie'i, 

so that the associated conformal map is 

z{t,w) = V2tw, \w\ = 1. 

On the other side, if we choose an initial contour to be different from the level set of 

the function A, the inverse problem may have a unique solution modulo an arbitrary 
function of A. Obviously, the simplest way to verify if the "movie" in the previous 
example is the Laplacian growth, is to run the "movie" again but with a non- 
circular initial configuration. Then the ambiguity of the previous example should 
disappear. But this way of removing non-uniqueness will require two "movies" . 
Thus, it is probably correct to say that there a continuum of A's, which correspond 
to the same "movie" , and that there are two sources of non-uniqueness: (i) any 
smooth function of A can replace the original function A(x) and (ii) any smooth 
function of p can be multiplied by the original A(x) without a change of the "movie" . 
This is an interesting fc^ature that merits further investigation. Here, we have simply 
pointed out (by constructing an example) the non-uniqueness of elliptic operators 
L providing the same evolution (i.e., a "movie"). 

7. SCHWARZ FUNCTION IN ELLIPTIC GROWTH 

It is instructive to state the moving boundary problem in terms of the Schwarz 
function S{z). Start with the equation for the velocity of the moving boundary 

V = diq. By Proposition 2.3 its left hand side equals St/{2iy/S^), while the right 
hand side can be written as de{C — p)/i = d^Q/i since dip = 0. (Here, as before 
C,=p + iq.) Then, 

Ct/i = {ztd.C + z,d,o/i = (d,c + s,d,o/iVs,. 

Hence, in virtue of the Beltrami equation (4.3) and the identity 

St 



we obtain the following evolution law for the Schwarz function. 
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Proposition 7.1. Under the above assumptions, 

(7.1) St{t, z) = 2(G + = 2d,az, S{t, z)) 
for all z in a neighborhood U ofr{t). 

Notice that the above formula is surprisingly similar to the evolution law of 
S{t, z) in the Laplacian growth process - cf. Theorem 2.5. 

7.1. Dynamics of singularities of the Schwarz function. The above proposi- 
tion has interesting consequences. We consider first the dynamics of the singulari- 
ties of the Schwarz function and associated universal quadrature formulas that are 
preserved during the elliptic growth. We have seen in the section devoted to the 
Laplacian growth that the poles, or more generally, the Cauchy integral density of 
the Schwarz functions of the moving boundaries are unchanged with the exception 
of the residue of the pole at z = which depends linearly on time t. Naturally, we 
expect that the relations in the elliptic growth process are more complicated. The 
present subsection collects some observations along these lines. 

Theorem 7.2. Let f{z) be a real analytic function defined in a neighborhood of 
the closed domains D(t) which are moving according to the elliptic growth law with 
an associated operator L. Then, 

(7.2) ^ / fdArea = 27r/(0), 

"I JD{t) 

where f solves the elliptic Dirichlet problem: 

Lu = in D{t), u\r(t) = /• 

Proof. Using the notations introduced in the earlier sections and previous compu- 
tations, we infer 

^ / fdArea = [ fVdi = [ f Xdnpdi = 2n f (0) . 

□ 

If the boundaries T{t) remain real analytic during the time of growth as is tacitly 
assumed throughout this note we can, as before, decompose the Schwarz function 
as follows: 

S{t,z)^S+it,z)-S-it,z), zer{t), 

where 

S.it,z) = -[ z^K. 
Jk o- - ^ 

Here, K is a compact subset of D{t) independent of t and ;U is a complex valued 
measure smoothly depending on t. Then, assuming in addition that / is (complex) 
analytic in a neighborhood of D{t), we find 
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This representation of the derivative of the average of an analytic function be- 
comes interesting in several particular cases. In particular, as in the conservation 
law obtained by taking / — 1, cf. Corollary 2.7 and Corollary 3.2, it follows that 

K dt 

Another application is given by the following. 

Proposition 7.3. Assume that for all times t the singularities of the Schwarz 
function contained in D(t) are simple poles; i.e., we have in D{t): 

N ^ 

S-{t, z) — — ^ irtS ^ (analytic remainder), 

with bj{t) G D{t) for all j and t. Then, 

N 

(7.3) 5^K.(t)/(6,(t)) + a,{tmt)f'{b,m = 2/(0) 

j=i 

for every analytic function f . 

Proof. According to the computations in the proof of Theorem 7.2, we have 



IJ,{t,a) 



d[J2f=iO'j{t)Sb^(t)]i<7) 



hence 



d r ^ 

- / fdArea - 7r^[a;.(i)/(6, + a,(i)6;.(t)/'(6,(i))] 



dt 

and the statement follows from (7.2). □ 



Now, by choosing 2A'^ linearly independent analytic functions fk we can, based 
on (|7.3p . form a linear system of equations that at least in principle determines the 
dynamics of the poles bj and residues Oj. 

8. Herglotz theorem and generating closed-form solutions 

The Herglotz theorem [JU] establishes a fascinating one-to-one correspondence 
between the singularities of the Schwarz function of the contour and the singularities 
of conformal maps from a vicinity of the unit circle to a vicinity of the contour under 
consideration. Basically, it states that if a is a singular point of a conformal map 
f{w) from the unit disk © to the domain D such that dD = {/(e''?);g G [0,27r]}, 
then the Schwarz function of the curve dD has the singularity of the same kind at 
the point 

(8.1) b=f{l/d). 

Moreover, if an isolated singularity (i.e., a pole, an algebraic singularity, or a loga- 
rithmic singularity) at a appears in the function f{w) with a coefficient A, then the 
corresponding singularity b is present in the Schwarz function with the coefficicnti? 
determined from 

(8.2) B^Ai-a^f'{l/a)r, 

where m is the multiplicity of a pole if a is a pole; is a rational number if a is a 
an algebraic branch point; or, if a is a logarithmic singularity, m is equal to zero. 
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Actually, the last two equations, which we call Herglotz' theorem, follow easily from 
the representation of the Schwarz function S{z) in terms of the conformal map / 

s{z)^fo{i/r'){z). 

There is some evidence that Herglotz' theorem should be helpful in solving the 
elliptic growth problem in terms of z = f{t, w) and generating exact solutions in 
the closed form. Here, we present a naive sketch of how we might expect generating 
of exact solutions for the elliptic growth should work. 

First, find ^ = p + as a function of z ^ x + iy and z — x — iy, either by solving 
the Beltrami equations (|4.3p or (|4.4p for a given A and a given initial domain D{0), 
or by solving the Dirichlet problem in -D(O), thus finding the Green function p and 
consequently calculating the conjugate function q from the generalized Cauchy- 
Riemann equations (4.1), (4.2). Then, the equation z) = implicitly defines 
the Schwarz function of the moving boundary. 

Second, either substitute ({z, S{z)) into the equation (|7.ip for the dynamics of 
the Schwarz function S{z) (scenario A), or, "if one gets lucky", try to invert the 
solution C(z,z), thus obtaining the function z(C,^) (scenario B). 

Scenario A: As the third step, identify singularities of S{z) (and their dynamics) 
that are already built in (|7.1|) through the function ({z, S{z)) found at the previous 
step and, also, identify the singularities of S{z) which are not the singularities of 
the RHS of (|7.1|) . The latter are time- independent as one can see from (|7.ip and 
represent constants of motion associated with the dynamics of the growth. 

Scenario A. Fourth step. Using the formulae (|8.ip . (|8.2p provided by the Herglotz 
theorem recover an explicit form for the moving boundary, z = /(i, e**), with the 
time dynamics of all parameters of / given implicitly by (|8.ip . (|8.2p . 

Scenario B. Third step. Restrict z(C, C,) to the imaginary axis of C, that is the 
axis p = 0, thus defining the function /(e**) = C,{iq,—iq), which, as was shown 
above, satisfies the equation p.2p for area preserving diffeomorphisms. 

Scenario B. Fourth step. Substitute /(e*'') into (|1.2p as an initial condition with 
time-dependent parameters and solve it. (The technique of integration is the same 
as that for the Laplacian growth [41], but the singularities of / are now lie both 
inside and outside the unit circle.) Alternatively, find the singularities of S{z) that 
correspond to the singularities of /(z) through (|8.ip . (|8.2p and calculate their time 
dynamics using (|7.ip . In either case the Herglotz theorem plays the central role in 
finding exact solutions for elliptic growth. 

While a more complete theory for solving (II. 2p for the elliptic growth exactly 
will be published elsewhere, we will present below two examples of exact solutions 
in the case when the singularities of f{w) are poles, both inside and outside the 
unit circle. 



8.1. Example: two simple moving poles. Let us take 
(8-3) ^(e^'') = re^' + + 

as an initial condition for the equation (|1.2p . assuming |ai| < 1 and \a2\ > 1. 
One can verify by a direct substitution that (|8.3p is a solution of p.2p with time 
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dependent poles ai, 02, residues Ai, A2 and the conformal radius r. The time 
dependence of these parameters is given by the equations 

(8.4) 5i,2 

(8.5) Bi^2 

(8.6) 2t + C 

where 61,2, -81^2, and C are constants of integration. Actually, the RHS of the last 
equation is the area of the domain D(t) enclosed by the contour z(i, e*^). When 
A2 = this is a standard Laplacian growth with a simple pole at ai (t) . In this case, 
the contour is known to develop a finite time singularity by forming a cusp. When 
Ai — 0, the process becomes a so-called inverse Laplacian growth, when a more 
viscous fluid displaces a less viscous one oppositely to a standard situation in which 
it is the other way around. The inverse Laplacian growth is stable and the shape 
rounds off during the evolution forming as a rule a circle as a long time limit. The 
case when both Ai and A2 are not equal to zero is a general case with a nontrivial 
dynamics caused by an interlay between the stabilization of a contour, due to the 
term in (|8.3p with the pole a2 that lies outside the unit circle, and destabilization, 
due to the term in (|8.3p with the pole oi inside the unit circle. 

It is interesting to note that the constants of integration in the LHS of (|8.4[) -f8.6) 
describe singularities of the Schwarz function of the moving contour, namely 61,2 
and Bi 2 are the simple poles and residues of S{z) respectively. Note that bi ^ D{t) 
while 62 G D{t). Finally, 2t + C is the residue of the Schwarz function at the simple 
pole at the origin, which represents the area of D(t). It only remains now to find 
the explicit expression of the function A(x) for this process described by (|8.3p . 

8.2. Example: two multiple stationary poles. For the second example let us 
take, as the initial condition for the (2.1), a function 

(8.7) zie'") = re'" + ae^^^"")' + 6e^(i+")«, 

that describes a contour with 71-fold symmetry which represents a circle of radius 
r modulated by a monochromatic wave in such a way that exactly n waves of an 
amplitude •y/2(|ap -I- |&p) fit the circumference (at least when \a\ and \b\ are both 
small with respect to r). 

One can substitute (|8.7p into (2.1) and verify that (18. 7p is a solution of the area 
preserving diffeomorphism (2.1) if the time dependent parameters a, 6, and r obey 
the following algebraic equations: 

(8.8) A = ai/(i-")6i/(i+"), 

(8.9) - = (1 + n)a^/^''-^^ + (1 - n)Aa^/^^~''\ 
r 

(8.10) 2t + C = - (n - 1) \af + {n + 1) \bf , 

where A, B, and C are constants of integration. 

The comments that can be made here about the dynamics of the contour de- 
scribed by (|8.7p are very similar to those made above for the first example. When 
6 = 0, the moving curve blows up in a finite time by forming cusps and thus ceases 
to exist after that. When a — the process is the inverse (stable) Laplacian growth, 
it smoothes the curve to a circle in a long term asymptotics. When both a ^ and 
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6^0, this is an intermediate case which can, for instance, stabihze the dynamics 
and prevent the finite time blow up with a proper choice of initial parameters. 

Just as in the previous example, the constants of integration describe singularities 
of the Schwarz function, namely A and B are the (1 — n)*'' and (l + n)*'* coefficients 
of the formal Laurent expansion of S(z); while 2t + C is the area of D{t) enclosed by 
z{t,e'"'). Unfortunately, again, it is not clear at the moment what elliptic parameter 
A is associated with this process. This question requires further investigation. 

9. Conclusions 

In this work we have shown that the elliptic growth processes, which present 
a natural generalization of the Laplacian growth, possess some remarkable math- 
ematical properties strongly resembling those pertinent to the Laplacian growth. 
Specifically, there is an infinite set of conservation laws; these conservation laws 
can be interpreted in terms of potential theory; there exists a parametrization of 
the moving interface by the stream function of an associated fluid velocity vector 
field and, last but not least, there are several interesting accompanying features 
related to the singularities of the Schwarz functions of the moving boundaries. We 
expect the latter to be especially helpful in generating and illuminating particular 
closed form solutions of elliptic growth problems. We think that the next step in 
this direction should be a search for a class of multipliers A, which will allow ex- 
plicit solutions of the Beltrami equations and, therefore, will offer new closed form 
solutions of the elliptic growth phenomenon. 
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